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Abstract 

Using a modified model of market dynamics with price expectations equivalent to a damped harmonic oscillator 
model, we analyze the reaction of financial markets to shocks. In order to do this, we gather data from indices of a 
variety of financial markets for the 1987 Black Monday, the Russian crisis of 1998, the crash after September 11th 
(2001), and the recent downturn of markets due to the subprime mortgage crisis in the USA (2008). Analyzing those 
data we were able to establish the amount by which each market felt the shocks, a dampening factor which expresses 
the capacity of a market of absorving a shock, and also a frequency related with volatility after the shock. The results 
gauge the efficiency of different markets in recovering from such shocks, and measure some level of dependance between 
them. 

1 Introduction 

Market crashes are a phenomenon that has been gathering much attention in the past four years, particularly due to 
the ongoing subprime mortgage crisis, but it is a subject that has been attracting the attention of numerous researchers, 
particularly economists, for the last few decades. The scientific literature shows many attempts at modeling the transmis- 
sion of volatility (contagion) between stock markets pQ- [5], how the correlation between them change with time [9"]-|ll|. 
and how it tends to increase in times of crisis [H]-[21]- Recently, physicists have been building models trying to explain 
how and why such crashes occur. One particularly successful model (although not completely), the log-periodic model 
[2S]-[55], predicts crashes as the possible outcomes of increasingly fast oscillations due to the interaction between agents 
in a market. 

In [53], we have shown how stock markets tend to behave similarly in times of crisis, responding to a strong world 
market movement that may be the result of news shared by all markets or back reactions to other market's movements. 
This market movement increases in times of crisis and one may model individual stock markets' movements as the response 
of those markets to an external force (the market). 

Our work uses a modification of a simple model of market dynamics with price expectations which can be put into 
direct analogy with a model of a damped harmonic oscillator subject to some external force in order to study the behavior 
of some stock markets chosen to represent different parts of the world. Periods of crashes make their study simpler 
by the fact that often the large oscillations of the markets overshadow their typical random walk behavior. Linking the 
oscillations after a crash to a model of a damped harmonic oscillator also makes it possible to analyze some of the market's 
properties, such as resistance to change, correlation to the other markets, and volatility. 

In order to do our work, we use data from stock market indices in 9 countries: USA, Hong Kong, Japan, Germany, 
UK, Brazil, Mexico, South Korea, and Australia, which, together, represent about 60% of the world financial market 
capitalization (around 28 billion dollars per month). Two indices from the New York Stock Exchange (Dow Jones and 
S&P 500) are used as benchmarks. The fits of data to the model are made using simple minimum squares techniques 
with an eye for the high nonlinearity of the function being calibrated. 



2 The choice of crises 

We now explain why we chose the crises of 1987, 1998, 2001, and 2008. Figure 1 shows the monthly average of the daily 
performances of the S&P 500 index of the New York Stock Exchange from the beginning of 1980 to the middle of 2010. 

*E-mail: leonidassj@insper.edu.br (corresponding author) 
t E-mail: italopf@al.insper.edu.br 



1 



The crashes of 1987 (Black Monday), 1998 (Russian crisis), 2001 (September 11), and 2008 (subprime mortgage crisis) 
are outlined in order to show the periods we are focusing in our study of how markets react to strong shocks. These 
crises have been chosen because they were severe, global, and because they represent a variety of types. The crisis of 1987 
was a surprise, probably caused by endogenous reasons; the 1998 Russian crisis was also abrupt, but followed by a quick 
recovery; the crisis of September 11 was an example of a crash that was caused almost completely by external news; and 
the crisis of 2008, which extended from 2007 up to 2009, was an example of a series of downfalls of the world markets. 
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Figure 1: S&P 500 index (monthly averages), from the beginning of 1980 to the middle of 2010, with the chosen crises highlighted. 

It is well known that financial markets behave in good approximations of efficient markets, and most of the theory 
in finance is based on this assumption. That may be true for normal periods, but the efficient market hypothesis turns 
out to be wrong when markets are experiencing crises. Looking at the closing indices of every day in which there was 
negotiation, we consider now the log-returns of the S&P 500 index, given by 



s t = H p t)-HPt-i) 



p t - p t -i 



(i) 



In order to best visualize the most extreme points, we next plot the log- frequency distribution (figure 2), defined as 

log-density = ln(l + density) , (2) 
where density is the frequency distribution of the log-returns, 
log-density 
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Figure 2: log-density of log-returns of the S&P 500 index, computed from the beginning of 1980 to the end of 2008. 

If the probability distribution function relative to the log-returns were a Gaussian, then we were not not expect large 
deviations from the average of the log-returns. Nevertheless, one can clearly see from figure 2 that there are strong 
deviations, mainly towards negative log-returns, that cannot be explained by a Gaussian distribution. In figure 2, a 
Gaussian distribution function should be represented by a parabola, what is clearly not the case, so financial crises should 
not be studied using the efficient market hypothesis. In figure 2, we circled the densitites relative to the four crises we 
are studying, showing that they correspond to diverse degrees of deviation from the norm. 
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3 The model 



Now we shall explain the kind of dynamic equations we are using in our model, and why, starting by connecting them 
with the equations of a well known model of economic dynamics. 

3.1 Market dynamics with price expectation 

We shall begin with a very simple model of price change due to market demands, which is used in textbooks on economic 
dynamics |54j . The model is based on the following two equations of demanded quantity Qd and supplied quantity Q s as 
functions of the price P of a commodity and its first and second order derivatives: 

Qd = a> - axP + a 2 P - a 3 P , (3) 
Q s = p + p x P - p 2 P + p 2 P , (4) 

where all constants are considered real and positive, and where the dynamics is given by 

P = KQd - Qs) , (5) 

where A > 0. Here P stands for the derivative of P with respect to time t. 

Equation §3§ expresses the assumption that demand falls when price rises, and that demand grows when the market 
perceives that prices have an upward tendency, and falls when this tendency is decreasing. In a similar way, equation 
(|4|) expresses that supply rises when price goes up, falls when there is a further tendency of the price rising, and goes up 
when this tendency is growing. The dynamic equation ([5]) says that, when demand is greater than supply, prices tend to 
rise, and when supply is higher than demand, prices tend to drop. 

Such a model is supported by some experimental evidence [5 5) and is considered appropriate for market dynamics, 
supposing there is not too much noise. It leads to the following differential equation for the price as a function of time: 

A(«! + fa)P + [1 - X(a 2 + p 2 )} P + \{a x + px)P = A(a - Po) ■ (6) 

This differential equation has three types of solutions, depending on the parameters of the model, specifically, on the 
combination 

A = 4A 2 (« 3 + /3 3 )( ttl + Pi) - [1 - X(a 2 + p 2 )f . (7) 



If A > 0, then 
being c\ and c 2 constants, 

and 



P(t) = ci e rit + c 2 e r2t + P* , (8) 

_ - [1 - \(a 2 + f3 2 )] -VA = - [1 - X(a 2 + p 2 )] + VA 

n 2X(a 3 +p 3 ) ' r2 2X(a 3 +p 3 ) ' [ > 

p* = 2°^° . (io) 

This solution describes exponencial growth or exponential growth or fall up or down to an asymptotic limit. 
The second solution occurs when A = 0: 

P{t) = c 1 e rt + c 2 te rt + P* , (11) 

where c\ and c 2 are constants, 



r _ -[1- X(a 2 +p 2 )} 
2X(a 3 + /3 3 ) 

and P* is given by (fTU|) . This solution broadly describes the same situation as the first type of solution. 
The third solution occurs when A < 0, and is given by 



(12) 



P(t) = ci e~ pt cos(wt - ip) + P* , (13) 
where c± (the amplitude) and ip (the phase) are constants, 

V^A _ Vl 1 " A («2 + P2)] 2 - 4A 2 (a 3 + /3 3 )(ai + Pi) 

W ~ 2(a 3 +p 3 ) ~ 2(a 3 +p 3 ) ' (15) 
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and P* is given by (jlOp . Depending on /?, this equation describes an oscillatory movement that diminishes exponentially 
with time or increases exponentially in time. 

The following three figures (figures 3a, 3b, and 3c) illustrate the three types of solutions. As can be seen, P* is the 
price of equilibrium for the three solutions, where Pq is the starting price. 




1 * t ' * t 1 * t 



Figure 3: a) example of solution for A > 0. b) example of solution for A = 0. c) example of solution for A < 0. 

So now one may ask why such a simple model should be introduced here, and what does it have to do with financial 
market crashes. If one looks at figures 3a and 3b, one can see that it resembles the behavior of market indices when there 
is a slow drop in prices, given or taken some noise. Now, the third figure ilustrates aproximately what happens to markets 
when there is a sharp drop (a crash), followed by oscillations due to speculation which further disappear into the usual 
noise close to a price P* that the market considers more appropriate. 

Our insterest in this article is to analyze the third type of behavior, and to gauge how each financial market behaved 
in previous crises when viewed through the light of a model similar to this, although the model we shall be using is not 
exactly this one, as shall be explained now. 

When looking at a differential equation like (JSJ, one must consider the homogeneous part of this equation (left side 
of that particular one) to describe endogenous effects, due to the nature of the markets. The nonhomogeneus part (right 
side of that equation) describes exogenous effects. What we shall do now is add another term, S e~ at , which describes an 
exogenous shock in the market that is very sharp in the begining, and that drops very rapidly with time. Such a shock 
tries to mimic the effect of some external information that may well be information about the drops of other financial 
markets or the effect of endogenous phase transitions, in a particular market. Doing so for differential equation (j6|), one 
obtains 



A(a a + py)P + [1 - X(a 2 + p 2 )]P + A(ai + fo)P = A(a - Pa) + 5 e~ at , (16) 

where all constants are considered to be real and positive. 

The solution for this differential equation when there is an oscillatory behavior is 

P(t) = a e~ pt cos(wt - ip) + P* + b e- at , (17) 

where 

b- 6 (1 

~ X(a 3 + f3 3 )a 2 - [I - X(a 2 + P 2 )} a + ay + Pi ' 

and the other coefficients are given by equations (ll"4"|) and flTSj) . 

Before using this model in order to gauge real financial data, it is best if one first sees it through the eyes of another 
theory that originally comes from physics. 



3.2 Damped harmonic oscillator 

The model we have just described can be put in direct analogy with another model of a damped harmonic oscillator that 
makes it easier to comprehend the solutions we have just shown in the last subsection. 

We start by considering a body suspended by a spring. This body is subject to two forces: the weight P and the 
traction T exerted by the spring. When both forces are evenly matched, the body stays at the equilibrium point x = 0. 
When an external force is applied to it, then the body will be displaced by a length x, and a restoring force F will appear. 

This type of force is described by Hook's law, 

F = —kx mx = —kx , (19) 

where a = x is the acceleration, m is the body's mass, and k is the spring's constant, which depends on how the spring 
deforms. The solution to this differencial equation is a combination of sines and cosines, which means it is a purely 
oscillatory solution. 
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If one adds the term —jx to the right side of equation (TT91 , then one obtains 



m'x — —kx — "fx , (20) 

which is the equation of a damped harmonic oscillator. The solutions to this equation range from exponencial decay to 
oscillations that decay exponentially, just like the solutions for the market model with price expectations. 
We now consider a model given by 

mP = -k(P- P*) -jP + Se- at , (21) 

where P is the price of a stock, P — P* is the distance between the price and a price P* that the market considers to be 
fair, and S e~ at is the shock we have introduced in the last subsection. This equation may be written like 

mP + ~/P + kP = kP* +5 e~ at , (22) 

which is essentially the same equation as (fTl)]) . and has the following possible solutions: 



P(t) = Cl e-* + c 2 e-* + P* +be-«\ n = Wl 2 ^ r 2 = rfHg , A > ; (23) 

2m 2m 

P(t)= Cl e- rt +c 2 te- rt + P* +be~ at , r = ^L , A = 0; (24) 

2m 



P(t)=c 1 e-f 3t cos{wt-tp)+P* +be~ at , = J- , w = ^ 4r " fc ^ , A<0, (25) 

2m 2m 

where ci, c 2 , and ip are constants and 

A = 7 2 - 4mfc , 6 = - . (26) 

ma 1 — 7a + k 

Solution (|25|) is particularly interesting to us. It describes a damped harmonic oscillation of the body attached to 
the spring when subject to an external force that is very strong at the first moment, and then diminishes quickly. The 
solution may be compared directly with solution (|17[) if one sets 

TO = A(a 3 +/? 3 ) , 7 = 1 - A(a 2 + /3 2 ) , k = \{a x + ft) . (27) 

Thus, one can now make analogies and have a more physical picture of the model of market dynamics. Coefficients 
«3 and (83, which are the influence the acceleration of prices have on the variation of price, may be seen as the "mass" of 
the system (seasoned by constant A). Mass has the effect of inertia: a system with high mass is more resistant to shocks, 
but is also slower to settle after it has been put into motion. Coefficients a 2 and /3 2 , which gauge the dependance of P 
on itself (when multiplied by A), are the equivalent to the dampening factor 7. The coefficients ol\ and /3i, which are 
related with the dependance of P with the price P, are essentially the constant of the spring, and determine how strong 
the spring is. 

In what follows, we shall consider the oscillatory solution to be given by 

P(t) = A + B e- at + C e~ pt cos{wt - ip) , (28) 

where A, B, a, C, (3, w and ip are all parameters determined by the experimental data. Constant A establishes the 
average axis around which the solution oscillates, and is also its assymptotic solution; constants B and a determine the 
strength of the initial shock and how fast it decreases; constant C determines the amplitudes of the oscillations; constant 
j3, how fast those oscillations diminish with time; constant w, the frequency, how fast it oscillates (related with volatility 
after a crash); and tp is the phase, which determines where are the peaks of the oscillations. The price P will be exchanged 
for indices of markets, like the Dow Jones (USA) or the Hang Seng (Hong Kong). 

This model has no less than seven parameters, which must be determined by real data on market indices, and most 
of them are linked in a nonlinear way. Section 5 explains how we attempted to do that. 



4 Relation with the log-periodic model 

Sornette and Johansen, followed by other collaborators, developed a model for the behavior of financial markets prior to 
and after crashes [3S] -[52], as well as other economic systems [SS], [57] • This model has been extensively tested among 
the years (for recent, yet unpublished articles, see [55]-[SI]) on a great number of markets, and has also received a solid 
theoretical standing based on the theory of phase transitions and finance theory, also receiving a certain number of 
criticism |62j . [63] . Their model, in its simplest form, can be represented by the following equations: 

P(t) = A + B(t c - tf + C(t c - tf cos (w ln(t c -t)-<f>) (29) 
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prior to a crisis and 



P(t) =A + B(t- t c f + C{t - t c f cos [w \n{t - t c ) - 0) 



(30) 



after a crash. In these equations, t c is the critical time at which the crash occurs, where the model goes to infinity and 
loses validity. Their model is measured in months and years, and the oscillations get stronger as time approaches the 
critical value. 

Using the change of variable 

r = ln(i c - t) , (31) 

into equation (|3T)]) . one has 

P(t) = A + Be/ T + Ce^ cos {wT~(j)) , (32) 

which can be readily compared with equation (1281) , with the important difference that we have two distinct values for the 
coefficients of the two exponencials. Considering now that their model measures time in terms of months or years and 
ours measures time in terms of a few days, one can use the expansion In a; ~ [x — k) + (1/2) (a; — k) 2 around a Xq — k to 
establish that 

t « t - t c , (33) 

so that our model can be seen as a kind of approximation of the log-periodic model for times close to t c . 

This is an important aspect of our model, and we shall go back to this comparison after we start analyzing our results 
coming from real data. 



5 Data 



Let's now explain how the data were treated. First of all, in order to compare the various indices, we worked with 
normalized values of them, given by the indices divided by their average along the period being studied. This was not a 
problem in our case, for we only considered very small time intervals. 

We next show how the fit to equation (|28|) was made. Since we have seven parameters to ajust, and the function to be 
ajusted is highly nonlinear, the straight use of the minimum squares method usually leads to disaster. The use of numercial 
methods is also not very reliable, since the landscape of the error function is rather complex and those algorithms tend 
to look for the closest minimum. We are going to use the index for the Bovespa (Sao Paulo stock exchange, in Brazil) 
during the crash of 1987 in order to illustrate the procedure we used, for the fit is relatively complex compared with some 
others. 

First of all, we tried to remove the dependance on A, what just corresponds to reescaling the data. Using least squares 
method, which corresponds to minimizing the error function 



i=i 



where P(i) is given by (f2"5)l and Pi are the real data, we first find 
A = 1 (as expected, since the data are normalized). Second, we 
removed a term B exp(-at), and minimized the difference of the real 
data with this term subtracted from it, resulting in B — 9.34 and 
a = 0.36. Third, we removed a term C exp(— fit) and the resulting 
term cos(wi — ip), obtaining C = 0.14, j3 — 0.22, w — 2.92, and 
ip = 1.52. 

From this starting point, we tried to adjust the parameters again 
independently, using least squares method, and then all the param- 
eters together, so that the theoretical curve resembles as best as we 
can the real data. The parameters for this curve are A = 0.92, 
B = 8.1, a = 0.67, C = 0.21, fi = 0.15, w = 1.14, and p = -9.1. 
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Figure 4: first fit of the Ibovespa for the crash of 1987. 



After that, the error function can be plotted in terms of each of the many parameters while keeping the remaining 
parameters constant, in the vicinity of the parameters we obtained. There are clear values for the best fits of parameters 
A, B, a, C, and (3, in the vicinity of the values that have been chosen. For the phase <p, the best fit oscillates, as was to 
be expected. Now, for the parameter w, the function is more complex, and we can see minima near 0.6 and 2.8, being 
2.8 a lower minimum. We shall not consider higher values for w, because for large enough values of the frequency, all 
parameters can be adjusted at will, but without economic meaning. 

After we fix the value of w, we may now do some fine tunning minimizing the error function for the remaining 
parameters. This implies in changes in the value of iu, which in turn leads to changes in the remaining parameters. This 
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process can go on for a while, and it normally tends to converge, as long as the variations of w are done separately. The 
two graphs bellow show the fits obtained when one leaves from the values w = 0.6 and w = 2.8. Leaving from w = 0.6, 
the best fit is A = 0.9373, B = 6.0237, a = 0.6976, C = 0.0541, j3 = 0.0414, to = 0.8016, and cp = -1.3423. Leaving from 
w = 2.8, the best fit is A = 0.9355, B = 6.2195, a = 0.6788, C = 0.1948, (i = 0.2308, to = 2.8555, and tp = 1.2748. 
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Figure 5: a) best fit for w = 0.8016. b) best fit for w = 2.8555. 

What one can see is that, for w = 0.8016, the best fit ignores the first oscillation and is better adjusted to the last 
ones. For w — 2.8555, the curve fits well the first oscillations, becoming increasingly innacurate for the last ones. 

Before going into the detailed analysis of each index, we must explain which indices we chose, and why we have done 
so. First, we chose two distinct indices of the New York Exchange, the Dow Jones and the S&P 500, in order to compare 
the same financial market when viewed by two different angles. It also serves for us to gauge our ability to correctly 
compute the parameters of the model. We also chose the Nasdaq, which reflects a more technologically oriented market. 
The Hang Seng of the Hong Kong Exchange, and the Nikkei, of the Tokyo Stock Exchange, were chosen in order to 
represent the Asian markets. We chose DAX (Germany) and FTSE (UK) as representatives of the European market, and 
Ibovespa (Brazil), and IPC (Mexico) as representatives of Latin America. Finally, we chose Kospi (South Korea) and 
ASX (Australia) as representing other emerging markets. 



5.1 The crash of 1987 - Black Monday 

Monday, October 19, 1987, was a day of shocked investors, baffled analysts, and despaired traders. The collapse started 
in Hong Kong, and swept its way towards the west. By the end of the month, markets throughout the world had fallen 
more than 30%, in average. The world took many years to recover. 

It is still not well known what were the reasons for the crash of October, 1987. Some point to a failure in the computer 
operated systems, some to panic generated by the fear of another 1929, some to external economic influences. What is 
known is that the efficient market hypothesis failed terribly in that day, and in the days that followed. 

Figure 6 shows the Dow Jones index from April/1987 to April/1988. The period from 10/13/1987 to 11/08/1987 is 
highlighted, showing the time span we are considering for our model. 
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Figure 6: S&P 500 index from April/1987 to April/1988. The period from 10/13/1987 to 11/08/1987 is highlighted. 

Figure 7 shows the graphs of the best fits that were obtained for the 1987 crash. The error function was calculated 
for the first 12 days (from 10/13/1987 to 11/08/1987), thus giving a maximum priority to the fitting of the first few 
oscillations after the crash. This is an example of an endogenous crash without a quick recovery: there were no news that 
triggered the downfall of all financial markets in the world, but just a sequence of small bad news, and the news of the 
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f) Dax (West Germany) - 1987. 



g) FTSE (UK) - 1987. 
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k) Asx (Australia) - 1987. 



The parameters are displayed together in the next table. 



Parameter 


Dow Jones 


S&P 500 


Nasdaq 


Hang Seng 


Nikkei 


Dax 


FTSE 


Ibovespa 


IPC 


Kospi 


Asx 


A 


0.97 


0.97 


0.94 


0.87 


0.97 


0.83 


0.91 


0.94 


0.35 


-0.02 


0.90 


B 


0.59 


0.82 


0.76 


4.09 


1.14 


0.58 


0.95 


7.01 


1.80 


1.07 


1.90 


a 


0.70 


0.68 


0.36 


0.51 


0.53 


0.11 


0.29 


0.70 


0.08 


0.00 


0.40 


C 


0.16 


0.20 


0.17 


0.56 


0.12 


0.07 


0.18 


0.33 


-0.18 


-0.02 


0.27 


P 


0.18 


0.17 


0.18 


0.26 


0.13 


0.08 


0.21 


0.33 


0.04 


0.01 


0.19 


w 


1.31 


1.23 


1.50 


1.78 


1.88 


1.61 


1.76 


3.10 


1.06 


1.04 


1.70 


<P 


3.04 


3.80 


5.23 


1.83 


2.29 


5.34 


0.92 


2.55 


6.17 


7.02 


6.20 



Table 1: parameters for the fits of the indices being considered (1987). 



One can see that the majority of parameters are quite similar, and that all markets, with the exception of South 
Korea, which seems not to have felt the crisis at all, behaved in similar ways during and immediately after the crash. 
The similarity of the coefficients for the Dow Jones and the S%P 500 indices, both from the NYSE, can be a guide to the 
precision of our fits. Parameter A seems significant up to the second decimal digit; parameter B is less accurate, with 
precision only up to the first decimal digit. The remaining paramaters all seem significant up to the second decimal digit. 
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The parameter A shows that the Hong Kong, West Germany, and Mexico suffered the most severe losses after the 
crash. Parameter B shows the strength of the initial shock to each market. Hong Kong and Brazil seem to have felt 
the strongest shocks. Parameter a shows how fast the shocks were absorbed. The initial shock fell less steeply for West 
Germany and the UK. 

Parameters B and /3 show how fast the initial oscillations declined, and the initial amplitude of these oscillations. 
They were strongest for Hong Kong, but also diminished faster in that market. 

One key parameter to be looked at is the frequency w of the markets. They show that all markets, with the exception 
of Brazil, reacted with the same volatility just after the crash. This is a tendency that endured through the crises we are 
studying. The phase tp just shows how long it took for each market to be struck by the severe after-crash oscillations. 
Since the South Korean market has such distinct parameters, we are not considering this market in our analysis. 



5.2 The crash of 1998 - Russian crisis 

Following the so called Asian crisis, in 1997, the prices of commodities fell worldwide. Russia, whose economy is heavily 
based on exports of commodities, was much affected by that. Other internal factors, like the war in Chechnya, drove the 
Russian economy to the brink of collapse. As many countries had money invested in Russia, the crisis spread to foreign 
markets and caused a crash from which they recovered slowly. 

Figure 8 shows the Dow Jones index from February/1998 to February/1999, with the period that goes from from 
08/25/1998 to 10/14/1998 highlighted, which is the time span we are considering for our model. 
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Figure 8: S&P 500 index from February/1998 to February/1999. The period from 08/25/1998 to 10/14/1998 is highlighted. 

The graphs of the best fits that were obtained for the 1998 crash, from 08/25/1998 to 10/14/1998, are shown now. 
Here, we have enforced the restriction that /3 > 0. Had we not done so, then many of the oscillations would grow out of 
control, for the second oscillations after the crash were usually larger than the first ones. 
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Figure 9: a) Dow Jones - NYSE (USA) - 1998. 
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k) Asx (Australia) - 1998. 



The parameters are displayed together in the next table. 



Parameter 


Dow Jones 


S&P 500 


Nasdaq 


Hang Seng 


Nikkei 


Dax 


FTSE 


Ibovespa 


IPC 


Kospi 


Asx 


A 


0.98 


0.97 


0.98 


0.97 


1.03 


1.06 


1.02 


1.03 


0.92 


1.00 


1.00 


B 


0.24 


0.24 


0.40 


0.25 


0.24 


0.29 


0.19 


0.38 


0.23 


0.00 


0.80 


a 


0.56 


0.59 


0.78 


4.23 


1.77 


0.83 


0.74 


0.84 


1.00 


0.68 


5.18 


C 


0.02 


0.03 


0.05 


0.06 


0.02 


0.02 


0.01 


0.03 


0.03 


0.01 


0.04 


13 


0.00 


0.03 


0.10 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.17 


w 


1.66 


1.65 


1.62 


0.72 


1.37 


1.50 


1.64 


2.68 


2.64 


1.34 


0.67 




4.77 


4.77 


4.65 


1.64 


2.63 


4.36 


3.61 


4.71 


4.25 


4.29 


0.96 



Table 2: parameters for the fits of the indices being considered (1998). 



The first striking feature is that most of the parameters (3 are set to zero. This was built into the fitting process, for 
the oscillations just after the crash actually grew stronger in most markets, and that would ruin the approximation for 
later periods. Most countries showed faster recoveries than in 1987, what is shown by the larger values of parameter A. 
Brazil and Mexico had larger volatilities after the shock, followed by large drops in their indices. In the case of Brazil, 
the crisis was followed by strong speculative attacks against the local currency, and the market was very insecure at the 
time. Australia seemed not affected by the international financial crisis. 



5.3 The crash of 2001 - September 11 

In the morning of September 11, 2001, the USA suffered the most severe terrorist attack in all its history. The death toll 
was close to 3,000, including the terrorists who perpetrated the attack. The world watched in horror as two airplanes were 
shown live colliding with the twin towers of the World Trade Center, in New York. On that day, panic striken traders 
tried to get their money out of the stock exchanges and place their money in safer investments. Trade was cancelled in 
both the NYSE and the Nadaq. Markets recovered fast after the shock, for the economy had not been hit significantly 
and the stock market was healthy. 

The S&P 500 index from March/2001 to March/2001 is shown in figure 10. The period from 09/10/2001 to 10/22/2001, 
which we shall analyze, is highlighted. 

* S&P 500 
1300- - 

/ 

1200-- 
1100- - 
1000- - 

900 \ 1 1 1 1 1 1 *■ month 

3/01 5/01 7/01 9/01 ll/Ol 01/02 03/02 

Figure 10: S&P 500 index from March/2001 to March/2001. The period from 09/10/2001 to 10/22/2001 is highlighted. 

We next show the graphs of the best fits that were obtained for the 2001 crash, from 09/10/2001 to 10/22/2001, 
together with the appropriate parameters. September 11 was a fine example of a crash triggered by some external news, 
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without reasons for it of its own. In what follows, note that nearly all of the indices had a fast recovery after the first 
drop and the first two oscillations. 
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k) Asx (Australia) - 2001. 
The parameters are displayed together in the next table. 



Parameter 


Dow Jones 


S&P 500 


Nasdaq 


Hang Seng 


Nikkei 


Dax 


FTSE 


Ibovespa 


IPC 


Kospi 


Asx 


A 


0.94 


0.94 


0.93 


0.95 


0.97 


0.93 


0.43 


0.98 


0.95 


0.50 


0.96 


B 


0.45 


0.32 


0.30 


0.18 


0.12 


6.44 


0.57 


0.92 


0.39 


0.56 


0.14 


a 


1.03 


0.99 


0.67 


0.64 


1.46 


1.27 


0.00 


1.21 


1.51 


0.01 


0.44 


C 


0.04 


0.02 


0.02 


-0.03 


0.03 


-0.05 


-0.03 


0.06 


0.02 


-0.01 


0.02 


P 


0.06 


0.01 


0.03 


0.00 


0.06 


0.00 


0.09 


0.08 


0.00 


-0.03 


0.06 


w 


1.36 


1.29 


1.60 


1.88 


1.98 


1.10 


2.24 


0.85 


1.41 


1.06 


1.80 





4.62 


4.26 


5.77 


4.04 


1.53 


3.09 


2.78 


-0.10 


4.41 


7.44 


-0.04 



Table 3: parameters for the fits of the indices being considered (2001). 



What we have here is a much faster dissipating shock (higher values of a with small oscillations in its aftermath and 
small values of B). The most volatile markets after the shock were Japan and the UK (higher values for w). The low 
volatility for Brazil is mainly due to a choice between two values of w, one more volatile and the other less volatile, in 
the fitting process. 



5.4 The crash of 2008 - subprime mortgage Crisis 

We now analyze a somewhat different type of crash which occured recently, triggered by a subprime mortgage crisis in the 
USA. This crash is different because it is much slower than the others, with a mixture of the recipe we have seen in the 
other crashes. In order to describe this difference, let us consider the Dow Jones index from 09/12/2008 to 10/27/2008 
(figure 12). The highlighted are shows a period where the index fell by a large amount, but one can easily notice the 
market was already going down before that. 



S&P 500 




month 

3/08 5/08 7/08 9/08 11/08 01/09 03/09 
Figure 12: S&P 500 index from March/2008 to March/2008. The period from 09/12/2008 to 10/27/2008 is highlighted. 

The next figure shows the graphs with the best fits for 2008, from 09/12/2008 to 10/27/2008. 
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d) Hang Seng (Hong Kong) - 2008. 
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h) Ibovespa (Brazil) - 2008. 



16 



V 





5 10 15 



k) Asx (Australia) - 2008. 
The parameters are displayed together in the next table. 



Parameter 


Dow Jones 


S&P 500 


Nasdaq 


Hang Seng 


Nikkei 


Dax 


FTSE 


Ibovespa 


IPC 


Kospi 


Asx 


A 


0.89 


0.88 


0.89 


0.17 


0.87 


0.90 


0.86 


0.87 


0.88 


-1.97 


0.90 


B 


0.50 


0.52 


0.41 


0.90 


0.46 


0.54 


0.28 


0.50 


0.26 


3.05 


0.25 


a 


0.52 


0.54 


0.47 


0.01 


0.42 


0.58 


0.21 


0.42 


0.21 


0.00 


0.20 


C 


0.06 


-0.06 


0.06 


-0.05 


0.06 


-0.14 


-0.06 


0.04 


0.03 


-0.02 


-0.03 


P 


0.01 


0.01 


0.01 


0.00 


0.01 


0.06 


0.00 


0.00 


0.00 


0.00 


0.00 


w 


1.13 


1.13 


1.24 


1.51 


0.93 


0.96 


1.52 


1.22 


1.45 


1.52 


1.17 




4.05 


0.87 


4.90 


-0.47 


3.79 


0.64 


6.07 


5.21 


8.21 


-1.60 


2.48 



Table 4: parameters for the fits of the indices being considered (2008). 



Once more, it was necessary to hx /3 as being positive only, what drove some of its values to zero. Note that the values 
of A are low, what signals a tendency of the market to continue its way down after this shock. The post-shock oscillations 
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are all small, as shown by the small values of the paramater C, and do not dissipate fast (parameter j3). The volatility 
now is nearly the same for all markets, what may be a hint that the markets being studied are now more mature. 



5.5 Multiple shocks versus secondary effects 

An analysis of the data prior to the first day studied for the crisis of 2008, now from 09/12/2008 to 10/27/2008, reveals 
(figure 14) that the best fit for this sample of data is a damped harmonic oscillator with longer wavelength, which means 
a lower value for the frequency w. 




A 


-7.8 


B 


9.01 


a 


0.00 


C 


0.06 





0.01 


w 


0.37 


<P 


4.01 


E 


0.02 



5 10 15 20 25 30 
Figure 14: Dow Jones - NYSE (USA) - 2008. Fit with long wavelength. 

This behavior is mimicked by all the other markets we are studying. Note, however, that this fit does not capture the 
minor oscillations within the wave described nor the larger fall that occurs at around day 14 (10/01/2008). 

In order to try to make our model more adequate, we shall introduce a second shock, which occurs at time to, and that 
makes the Dow Jones index drop with greater intensity. Such a shock can be written as ee~ ( ' t H(t — to), where H{t — to) 
is the Heaviside function, given by 

*(*-*>) = {!!;!>£; ^ 

so that the shock only begins at time to, which shall be considered a new parameter of the model, along with parameters 
e and £• 

Now, from the point of view of damped harmonic oscillators, model (f2"2")) becomes 



mP + jP + kP = kP* + S e- at + e e- ct H(t - to) . 
This differential equation can be solved using the Laplace transform. The oscillatory solution is 



P{t) = A + Be~ at 
where 



Ce 



s(wt -ip) + D e - c(t - to) H(t -to) + E e -^(*-*o) cos [w(t -to)-r)) H (t - t a ) 



A = P* , B 
tp = arctg 
D 



ma 2 — 7a + k 
s + l + ( 7 /m)(P -P*) 
B 



7 

2m ' 
+ a-j3 



\J Amk — 7 2 



2m 



C 



B 



COS if 



mC, 2 — 7C + k 



V = — arctg (£ — 0) , E 



D 

cosry 



(36) 
(37) 

(38) 
(39) 
(40) 



Fitting the new parameters, as well as refitting the previous ones, we obtain the following curve. 



18 



V 



1.2 ■ 



1.1 



1 -- 



0.9 -- 



A 


0.49 


B 


0.70 


a 


0.02 


C 


-0.10 


P 


0.03 


w 


1.12 


f 


3.01 



D 


0.32 


E 


0.03 


c 


1.15 


V 


3.95 


to 


13.55 


E 


0.08 



0.8- 



-+- 



-+- 



-+- 



t 



5 10 15 20 25 30 
Figure 15: Dow Jones - NYSE (USA) - 2008. Fit with two shocks. 

One can see the fit is not very satisfying. We shall try, then, a different approach, using a second mode of vibration 
and trying to fit it on the previous data. 
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Figure 16: Dow Jones - NYSE (USA) 
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Fit with two modes of vibration. 
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B 
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0.02 


a 


0.00 




p2 


0.00 


c x 


0.06 




W2 
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Pi 
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<P2 
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Wi 
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E 
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The new model may be described by the following function: 

P{t) = A + B e- at + Ci e- pit cos{ Wl t -tpi) + C 2 e~ ht cos{w 2 t - tp 2 ) ■ (41) 

Although it has less parameters than the multiple shocks model, it seems to describe best the behavior of this particular 
market. Such function is typical of coupled harmonic oscillators, resulting from the interaction of one or more markets 
and subject to a strong, rapidly dissipating shock. So, what this suggests is that the crisis of 2008 was most likely not the 
result of the response of individual markets to common, multiple external shocks, but might be the result of backreactions 
among markets. 

We next make some comments on how financial markets are integrating in the past decades. This discussion is useful 
in our future attempts at linking the markets using an improvement of our model that acomodates many markets arc 
coupled, damped harmonic oscillators. 



6 Correlation 

One good measure of how different markets are correlated during a period in time is the average of their correlation 
matrix. In the crisis of 1987, the average correlation between the indices being considered here during the period of crisis 
was < C >= 0.32. 

Another measure of how correlated different markets are is obtained by calculating the eigenvalues of the correlation 
matrix between the indices of those markets. The largest eigenvector corresponds to a "market mode", which expresses the 
common variation of the indices. In the case of the indices being considered for the crisis of 1987, the largest eigenvector 
of the correlation matrix between the many indices was A = 4.89. The graphs in figure 9 show that some of the indices 
reacted much more strongly than others. 

In the crisis of 1998, the average correlation between the indices being considered here during the period of crisis was 

< C >— 0.36, and the largest eigenvector of the correlation matrix between the many indices was A = 4.40. The graphs 
in figure 11 show that few of the indices reacted much more strongly than others. 

In the crisis of 2001, the average correlation between the indices being considered here during the period of crisis was 

< C >— 0.30, and the largest eigenvector of the correlation matrix between the many indices was A = 4.78. Figure 13 
shows good correlation between many of the indices. 
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In the crisis of 2008, the average correlation between the indices being considered here during the period of crisis was 
< C >= 0.46, and the largest eigenvector of the correlation matrix between the many indices was A = 7.71. Figure 15 
shows that all the indices considered reacted with the same amplitude, showing a strong dependance of each with the 
others. 

7 Mapping back 

We now try to map the coefficients we obtained back to the original coefficients, first of the damped harmonic oscillator 
model subject to an exponentially decreasing force, and then to the model of market dynamics with price expectations. 
Using equations (|2"5j) and (j2"6")l . we obtain the following relations between coefficients: 

P* = A , ci = C , 7 = 2m(3 , S = B(ma 2 - + k) , k = — + raw 2 , (42) 

4ro 

where P* would be the final value of the index if the market ever stabilized, C\ may be related with the index at time 
zero by the equations 

= P * - p ° ~ b Pq = p (0 ) , (43) 
cos(— <p) 

where S is the initial strength of the shock, and k is the constant of the spring, which sets how stiff the spring is. One 
can readily see that most equations depend on the same coefficient, m, which represents the mass of the system. 

This leads us to a problem, that we dubbed "the problem of mass": what should represent the "mass" of a financial 
market? A massive market would be harder to be put into motion, but also slower to stop. It would also oscillate less 
rapidly than a market with little mass, and could be a major influence if markets would be linked in a network. 

One of the possible answers would be that the mass of a market should be the amount of money it trades. We shall 
try this approach now, using as an example the data obtained for 2008. The following table contains the capitalization 
(in billions of dollars), of the 10 financial markets we are considering in our analysis at September, 2008, which is the 
time at which we do our analysis of that crisis. 



Market 


NYSE 


Nasdaq 


HK 


Jap 


Gcr 


UK 


Bra 


Mex 


SK 


Aus 


Capitalization 


13,046 


2,904 


1,614 


3,334 


1,352 


2,565 


888 


329 


656 


926 



Table 5: market capitalizations, in billions of dollars, at September, 2008. 



Considering the NYSE as having mass one, and writing the masses of other markets in terms of it, we then obtain the 
following coefficients. 



Market 


DJ 


S&P 


Nasdaq 


HK 


Jap 


Ger 


UK 


Bra 


Mex 


SK 


Aus 


m 


1.000 


1.000 


0.223 


0.124 


0.256 


0.104 


0.197 


0.068 


0.025 


0.050 


0.071 


p* 


0.944 


0.945 


0.930 


0.947 


0.974 


0.927 


0.431 


0.983 


0.955 


0.503 


0.893 


Cl 


0.039 


0.025 


0.025 


-0.027 


0.033 


-0.053 


-0.033 


0.058 


0.024 


-0.010 


0.273 


7 


0.125 


0.026 


0.014 


0.000 


0.030 


0.000 


0.036 


0.011 


0.000 


-0.003 


0.023 


k 


1.856 


1.671 


0.567 


0.435 


1.003 


0.126 


0.992 


0.050 


0.050 


0.057 


0.206 


S 


1.254 


0.827 


0.201 


0.088 


0.181 


1.892 


0.567 


0.126 


0.042 


0.031 


0.389 



Table 6: coefficients of the damped harmonic oscillator model according to index (2008). Mass is chosen as capitalization. 

The results are not very clarifying, and this may result from the fact that the mass of the NYSE is much larger than 
most of the others, what leads to very small values for the dampening factor 7 and to the strength 5 of the initial shocks 
to some markets. 

Another approach to the problem of mass is to consider the inertial aspect of mass, which sets it as the difficulty in 
putting the market into motion. A variable that describes approximately how fast a market oscillates is its volatility, 
given by the absolute value of the standard deviation of its log-returns. Using the inverse of this parameter as mass, so 
that the mass of a market will be larger when volatility is small, and smaller when volatility is high, and setting this 
parameter as 1 for the NYSE during the period we are studying, we obtain the following table. 



Market 


DJ 


S&P 


Nasdaq 


HK 


Jap 


Ger 


UK 


Bra 


Mex 


SK 


Aus 


ra 


1.000 


0.915 


0.940 


0.837 


0.826 


1.081 


0.958 


0.686 


1.056 


1.023 


1.204 


p* 


0.944 


0.945 


0.930 


0.947 


0.974 


0.927 


0.431 


0.983 


0.955 


0.503 


0.893 


Cl 


0.039 


0.025 


0.025 


-0.027 


0.033 


-0.053 


-0.033 


0.058 


0.024 


-0.010 


0.273 


7 


0.125 


0.024 


0.058 


0.000 


0.097 


0.000 


0.175 


0.109 


0.000 


-0.063 


0.392 


k 


1.856 


1.530 


2.397 


2.944 


3.242 


1.318 


4.836 


0.499 


2.113 


1.150 


3.501 


S 


1.254 


0.757 


0.849 


0.598 


0.586 


19.734 


2.762 


1.271 


1.743 


0.640 


6.591 



20 



Table 7: coefficients of the damped harmonic oscillator model according to index (2008). Mass is chosen as the inverse of the 
average volatility. 

The markets now have masses that are closer to one another, and the results are more similar. The problema with 
this definition of mass is that it leads to very different masses in 1987, though, influencing the remaining coefficients. 

Another problem we have is that both definitions of mass change with time, and depend drastically on wether the 
stock markets are facing a crisis or not. We need a more stable definition for mass, and that is a topic for future research. 
That definition may come from the study of markets as a network of damped harmonic oscillators, and we hope we are 
able to reach it soon. 

8 Conclusion and future research 

Our work showed how stock markets may be modeled in periods of crisis as damped harmonic oscillators subject to an 
intense, but fast decreasing external force. The model may be put in analogy to a log-oscillator model for short time 
spans. The results show that some markets share some characteristics and differ in others like volatility after a crash, 
and absortion of the shocks generated by the crisis. It showed an increasing interaction between markets in times of crisis 
during the years via the correlation matrix of their log-returns. It also led us to believe that models of interacting markets 
as coupled harmonic oscillators may be a better way to understand the co-movements of financial markets in times of 
crisis. Our research shall concentrate on this topic from now on. 
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